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Abstract. Existence of the /f-dibaryon in the flavor SU (3) symmetric limit is studied by full QCD simulations on the lattice, 
in the approach recently developed for the baryon-baryon (BB) interactions. Potential of the flavor-singlet BB channel is 
derived from the Nambu-Bethe-Salpeter wave function, and a bound /f-dibaryon is discovered from it, with the binding 
energy of 20-50 MeV for the pseudo-scalar meson mass of 469-1 171 MeV. 
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The //-dibaryon, predicted by R. L. Jaffe in 34 years ago [1]. is one of most famous candidates of exotic-hadron. This 
prediction was based on the observations that the quark exclusion can be completely avoided due to the essentially 
flavor-singlet(M«c/(i55) nature of //-dibaryon, together with the large attraction from one-gluon-exchange interaction 
between quarks suggested in the quark model [1, 2]. Search for the //-dibaryon is one of the most challenging 
theoretical and experimental problems in the physics of strong interaction and quantum chromodynamics (QCD). 
Still it is not clear whether there exists the //-dibaryon in nature. Although deeply bound //-dibaryon with the binding 
energy Bh > 7 MeV from the AA threshold has been ruled out by the discovery of the double A hypernuclei, ^He [3], 
there still remains a possibility of a shallow bound state or a resonance in this channel [4]. 

Several lattice QCD calculations on H have been reported as reviewed in [5] (see also recent works [6, 7, 8]). 
However, there is a serious problem in studying dibaryons on the lattice: To accommodate two baryons inside the 
lattice volume, the spatial lattice size L should be large enough. Once L becomes large, however, energy levels of 
two baryons become dense, so that isolation of the ground state from the excited states is very difficult (quite a large 
imaginary-time t is required). All the previous works on dibaryons more or less face this issue. In this paper, we 
employ our original approach distinct from one in previous studies, and search the //-dibaryon in lattice QCD [9]. 
Because the flavor SU (3) breaking complicate calculation, we start with the flavor SU (3) limit. 



To study interaction or property of a hadron systems in lattice QCD, we utilize the Nambu-Bethe-Salpeter (NBS) wave 
function and the corresponding potential, instead of the energy eigenstate. Within the non-relativistic approximation, 
potential of a two-body interaction can be defined from the NBS wave function ^(r, t) through the Schrodinger 
equation in the Euclidean space-time [10, 11, 12], 



where a non-local but energy-independent potential U(r,f) is introduced, and jx is the reduced mass of the system 
with masses M\ and M^. The potential U(r,r l ) can be determined by solving this equation by using data of ^(r,?) 
measured on the lattice. Since only a single data of ^(r, t ) is available from one simulation, we cannot construct whole 
the potential U (r, r'), but can obtain the leading term of the velocity expansion of it. We describe the potential V (r). As 
shown in ref. [9], one obtains a sink-time independent V (r) even without isolating energy eigenstates. It is also shown 
that we can obtain a volume independent V{r), if we setup the volume larger than the minimum to accommodate the 
interaction, e.g. L c± 4 fm for the nucleon-nucleon interaction at M n > 450 MeV. Once we obtain a volume independent 
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TABLE 1. Summary of lattice parameters and hadron masses. 



a [fm] 


L [fm] 


Kuds 


Mp. s . [MeV] 


M vec [MeV] 


M hm [MeV] 


Af cfg 






0.13660 


1170.9(7) 


1510.4(0.9) 


2274(2) 


420 






0.13710 


1015.2(6) 


1360.6(1.1) 


2031(2) 


360 


0.121(2) 


3.87 


0.13760 


836.5(5) 


1188.9(0.9) 


1749(1) 


480 






0.13800 


672.3(6) 


1027.6(1.0) 


1484(2) 


360 






0.13840 


468.9(8) 


830.6(1.5) 


1163(2) 


600 



potential, any observable of the system, such as binding energy and scattering phase shift, can be obtained by solving 
the Schrodinger equation in the infinite volume. Note that, in contrast to the conventional Luscher's method [13, 14], 
we do not calculate the energy shift of two hadrons at finite L to access the observables in the infinite volume. Finally, 
the NBS wave function can be obtained by summing up the hadron four-pint function G4 which can be evaluated 
numerically in lattice QCD. 

LATTICE QCD SETUP 

For dynamical lattice QCD simulations, we have generated five ensembles of gauge configuration on 32 3 x 32 lattice. 
We have employed the renormalization group improved Iwasaki gauge action with the coupling constant j3 = 1 .83, and 
the non-perturbatively 0(a) improved Wilson quark action. Quark propagators are calculated for the quark wall source 
with the Dirichlet boundary condition in the temporal direction. The point type octet-baryon field operators are used 
at sink. The sink BB operator is projected to the A\ representation of the cubic group, so that the NBS wave function 
is dominated by the S-wave component. For the time derivative of the NBS wave function, we adopt the symmetric 
difference on the lattice. To enhance signal over noise, we average 16 measurements for each configuration, together 
with the average between forward and backward propagation in time. Statistical errors are evaluated by the Jackknife 
method. Lattice parameters such as lattice spacing a, lattice size L, the hopping parameter K u ^ s , the number of gauge 
configurations N C f g , together with the hadron masses are summarized in Table 1. 

RESULTS 

The left panel of Fig.l shows the flavor-singlet BB potential at each quark mass, extracted from our numerical 
simulation of full QCD. One see that it has an "attractive core" and its range is well localized in space. This entirely 
attractive potential prove that the quark model prediction to the flavor-singlet BB channel is essentially correct. One 
can see that the long range part of the attraction tends to increase as the quark mass decreases. 

By solving the Schrodinger equation involving this potential, we can obtain any information of the system, for 
instance, energy of the ground state. It turns out that there is only one bound state in each quark mass [9]. The 
right panel of Fig.l shows the energy and the root-mean-square (rms) distance of the bound state. The result for 
Mp.s. = 469 MeV is only preliminary because of insufficient number of independent gauge configuration. These bound 
state correspond to the //-dibaryon predicted by Jaffe. Our lattice QCD calculation shows that a stable //-dibaryon 
certainly exists in the flavor St/ (3) limit world with the binding energy of 20-50 MeV for the present quark masses. 
Despite that the attractive potential become stronger as quark mass decreases, the resultant binding energies of the 
//-dibaryon decrees in the present range of the quark masses. This is due to the fact that the increase of the attraction 
toward the lighter quark mass is compensated by the increase of the kinetic energy for the lighter baryon mass. The 
rms distance \J (r 2 ) is a measure of spacial distribution of baryonic matter in the //-dibaryon, which can be compared 
to the point matter rms distance of the deuteron in nature, 1.9 x 2 = 3.8 fm. One may be able to get feeling of the 
//-dibaryon from this comparison: //-dibaryon is much compact compared to usual nuclei. 

When we solve the Schrodinger equation, we use the potential expressed in terms of an analytic function fitted to 
the data. Error from the choice for the analytic function to fit is negligible (less than few %). Small systematic error 
arise from choice for the sink-time-slice. The final result of the binding energy of the //-dibaryon in the SU(3) limit 
described as §h becomes, for example, 



Mp. s . = 837 MeV: Bh = 37.8(3.1)(4.2) MeV 
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FIGURE 1. Left: Potential of the flavor-singlet BB channel. Right: The ground state of the flavor-singlet BB channel. 



with statistical error (first) and systematic error (second). 

One may think such a deeply bound //-dibaryon is ridiculous or ruled out by the discovery of the double A nuclei. 
The present Bh should be interpreted as the binding energy from the BB threshold averaged in the (S,I) — (—2,0) 
sector. The three BB thresholds in this sector split largely in the real world with the flavor 51/(3) breaking. Therefore, 
we expect that the binding energy of the //-dibaryon in the real world, measured from the AA threshold, is much 
smaller than the present value or even //-dibaryon goes to above the AA threshold. To make a definite conclusion on 
this point, however, we need (2+l)-flavor lattice QCD simulations and the analysis in AA — NZ — EE coupled channel. 
Study along this direction is in progress [15]. 
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